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Abstract. The aim of this paper is to test numerically the predictions of the Mode 
Coupling Theory (MCT) of the glass transition and study its finite size scaling 
properties in a model with an exact MCT transition, which we choose to be the 
fully connected Random Orthogonal Model. Surprisingly, some predictions are verified 
while others seem clearly violated, with inconsistent values of some MCT exponents. 
We show that this is due to strong pre-asymptotic effects that disappear only in a 
surprisingly narrow region around the critical point. Our study of Finite Size Scaling 
(FSS) show that standard theory valid for pure systems fails because of strong sample 
to sample fluctuations. We propose a modified form of FSS that accounts well for 
our results. En passant, we also give new theoretical insights about FSS in disordered 
systems above their upper critical dimension. Our conclusion is that the quantitative 
predictions of MCT are exceedingly difficult to test even for models for which MCT is 
exact. Our results highlight that some predictions are more robust than others. This 
could provide useful guidance when dealing with experimental data. 



PACS numbers: 64.70.Q,64.60.Ht,64.60.F 
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1. Introduction 

In spite of all its shortcomings, the Mode-Coupling Theory (MCT) [TJ [2] of the glass 
transition has most valuably contributed to our understanding of the slowing down 
of super-cooled liquids. It was the first ab initio theory to make precise, quantitative 
predictions about the appearance of a two-step relaxation process in super-cooled liquids 
and hard sphere colloids as these systems become cooler or denser. MCT predicts in 
particular a non trivial "/3-relaxation" regime where dynamical correlation functions 
pause around a plateau value before finally relaxing to zero. Around this plateau value, 
power-law regimes in time are anticipated, together with the divergence of two distinct 
relaxation times, r a and Tp, as the MCT critical temperature Td is approached. Although 
this divergence is smeared out by activated events in real liquids, the two-step relaxation 
picture suggested by MCT seems to account quite well for the first few decades of the 
increase of r a [H 12] ■ 

Quite remarkably, it was observed soon after by Kirkpatrick, Thirumalai and 
Wolynes J3j Hj, that the MCT equations in fact describe the exact evolution of the 
correlation function of mean-field p-spin and Potts glasses [HJ E]- Furthermore, some 
physical properties of these p-spin glasses (such as the existence of an entropy crisis at 
a well defined temperature T s < T d ) are tantalizingly similar to those of super-cooled 
liquids, even beyond the MCT regime. This has led to: 

• A better understanding of the physical nature of the MCT transition, in terms 
of a "demixing" of the unstable saddle points of the energy landscape and the 
(meta) stable minima: above a certain energy threshold E t h only the former type is 
found, while minima all sit below this threshold (for a recent review, see [7]). 

• The elaboration of a consistent phenomenological description of the glass transition, 
called the "Random First Order Transition" (or RFOT) theory by Wolynes and 
collaborators [3j HJ E] , within which the activated processes allowing the system to 
hop between metastable minima is given a precise interpretation in terms of spatial 
rearrangements. 

More recently, it has been argued that MCT can be thought of as a Landau theory of 
the glass transition, where the order parameter is the (small) time dependent difference 
between the correlation function and its plateau value [9]. It was furthermore shown 
that the MCT transition is accompanied by the divergence of a dynamical correlation 
length [TOj, [TTJ (see also [12]), which gives a quantitative meaning (within MCT) to the 
concept of heterogeneous dynamics. Correspondingly, critical fluctuations are expected 
close to the MCT transition Td, and become dominant in spatial dimensions d < d u , 
where the value of the upper critical dimension is d u = 6 or 8 depending on the existence 
of conserved variables (see [131 EH])- For the model considered below, d u = 6. 

In view of the central role of the Mode-Coupling Theory in our current 
understanding of the glass transition, it is somewhat surprising that so little numerical 
work has been devoted to models for which the MCT equations are believed to be exact. 
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To our knowledge there is in fact no exhaustive treatment of a spin model displaying 
an exact MCT-transition. Such a study is valuable for several reasons. First, it is 
important to know how precisely the MCT predictions can be tested on a model where 
the theory is supposed to be exact, and for which all the excuses for MCT's failures in 
real systems (uncontrolled approximations, activated events, low dimensions, etc.) are 
absent. Second, the detailed study of finite size effects is important, since one expects in 
that case to observe in a controlled way the famous crossover between the MCT regime 
and the activated regime. Furthermore, since the short-range nature of the interactions 
in liquids should somehow lead to finite size corrections to the MCT equations, the 
results of this analysis should bring important insights, and perhaps help to understand 
the somewhat unexpected results of Karmakar et al. [15J. 

One reason explaining why these numerical studies are scarce is the slow dynamics 
of these models. Even more so than in other spin glass models a good sampling requires 
a number of Monte Carlo iterations that grows rapidly with the system size even with 
an extremely efficient sampling algorithm, like the parallel tempering algorithm [16J. 
An other reason is that the paradigmatic p-spin model is a p-body interaction model, 
which only leads to MCT-like dynamics for p > 3. Altogether this leads to quite 
heavy simulations (see [E]). Other mean field models in the MCT class exist, like 
the fully-connected q-states Potts model (for q > 4) [3] or the random orthogonal 
model (ROM) [HI fl9] . which are two-body spin models. But a recent detailed study 
of the 10-states Potts model [201 ED I22J [23] has led to rather strange results: absence 
of the "cage effect", unusual finite-size scaling exponent. Nevertheless some behavior 
reminiscent of the MCT-transition seems to emerge. The ROM case is even worse; 
there is no precise study in the temperature range close to the dynamic temperature, 
where all previous simulations have fallen out-of-equilibrium even for rather small system 

(N = 50 pa I2H 1251 ESS]). 

Our project is to provide an exhaustive numerical analysis of the statics and 
the dynamics of a model for which the Mode Coupling Theory is exact. We want 
in particular to : (i) test numerically the MCT predictions concerning the two-point 
relaxation function (<&(£)) and the four-point correlation functions (Xi(t)) th & t describe 
dynamical heterogeneities, (ii) study the pre-asymptotic corrections to the critical 
behavior and (iii) analyze the finite size scaling of the MCT transition. 
To achieve such a program, one needs a two-body model with well-separated static 
(T s ) and dynamic (To) transition temperatures. We have found that a certain variant 
of the ROM satisfies these constraints. We also need an efficient algorithm, since the 
ROM turns out to be an extremely difficult model to simulate. If one uses a simple 
algorithm like thermal annealing, the convergence is so poor that the average energy 
never goes below E t h below T^, even for small system sizes. In this work we will use 
the best algorithm to date for spin glass numerical simulations, namely the parallel 
tempering algorithm [16]. In spite of tremendous numerical efforts, we are still limited 
to rather small systems (up to 256 spins). However, our simulations allow us to reach an 
interesting but somewhat unexpected conclusion: preasymptotic corrections and finite 
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size effects are so strong that a direct observation of the MCT predictions is extremely 
difficult. Great care must be exercised to extract meaningful value of critical exponents, 
even for a model for which MCT is in principle exact! Such problems are expected to 
arise in the analysis of experimental data as well, and we believe that our work highlights 
some caveats concerning the validity and relevance of MCT in practice. 

Our first main result consists in establishing that pre-asymptotic effects are 
extremely important until one approaches the MCT transition temperature up to 
fractions of percent! We show this effect in two ways: (a) our numerical simulations of 
the ROM unveil that several of the MCT predictions are still not verified quantitatively 
for temperatures a few percent higher than the critical one; (b) we solve numerically 
the schematic Mode Coupling Theory equations and show that very similar conclusions 
are also reached in this framework. Similar findings were obtained in J27J by comparing 
MCT to experimental data. 

Our second main result is a detailed study of finite size scaling of the MCT 
transition. This is both of practical and theoretical interest when one attempts to 
determine some MCT exponents. We find that sample to sample fluctuations play a 
crucial role. As we shall explain, thermal fluctuations would lead naturally to finite 
size effects that becomes relevant for distance from the critical temperature of the order 
iV~ 2 / 3 where N denotes the system size. However, disorder fluctuations, that is sample 
to sample fluctuations of the dynamic temperature, are of order N~ l l 2 . The latter 
fluctuations therefore dominate and lead to a very subtle finite size scaling behavior. 
Our study lead us to a generalization of the Harris criterion suitable for disordered 
systems above their upper critical dimension. 

The paper is organized as follows. After recalling the main predictions of the replica 
method (for the statics) and of the MCT (for the dynamics) in section |2j we discuss 
our numerical results for the statics of the ROM in section [3} In |4j we investigate 
the equilibrium dynamics of the model and compare the results with the predictions 
of MCT. We show that the glaring discrepancies come from the unexpectedly narrow 
critical window around the MCT transition point, which we cannot access numerically 
without large finite size effects. We then turn in section [5] to a detailed study of these 
finite size effects, and to the possibility of using dynamical finite-size scaling for this 
model. We find that naive finite-size scaling theory fails to account for our results, 
because of the strong sample to sample dependence of the critical temperature. We show 
how to understand phenomenologically our results in the rest of section [5j relegating to 
appendices more precise statements on the apparent breakdown of the Harris criterion 
in high dimensions, and the exact solution of the fully-connected disordered Blume- 
Capel model, which provides an explicit illustration of our arguments. The last section 
is dedicated to the conclusion and open questions. 
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2. 1-RSB models: A short summary of known results 

As mentioned in the introduction, the starting point of the Random First Order Theory 
of the glass transition is the analysis of the so-called discontinuous spin glasses, namely 
mean field disordered models with an order parameter that has a jump at the transition. 
In terms of the replica method, these are the ones solved by the so-called one-step replica 
symmetry breaking ansatz (1-RSB, see [28l [29] and references therein). Examples are 
the p-spin (spherical) model and the Random Orthogonal Model that we are going to 
study in detail in the following. 

The physical behavior of these models is particularly transparent in terms of the 
Thouless-Anderson-Palmer (TAP) (30] approach that allows to analyze the free energy 
landscape of the model. Technically, this is a Legendre transform of the free energy as a 
function of all the local magnetizations. The minima of the TAP free energy correspond 
to the thermodynamic states of the system, like the two minima of the Curie- Weiss 
free energy represent the two low temperature ferromagnetic states in the (completely 
connected) Ising model. 

For a 1-RSB system, the analysis of the number M of solutions of the TAP 
equations that contribute to the free energy density shows that there exists two transition 
temperatures: one static (T s ) and the second dynamic (T rf ), with T d > T s . When T > T d 
or T < T s , the complexity, defined as iV -1 In Af, vanishes in the large N limit (N is the 
total number of spins). Above T d , this is because the energy is higher than the threshold 
value E th mentioned in the introduction, so that typical states are unstable saddles, 
hence the system is in a paramagnetic state. Below T s , this is because the low lying 
amorphous minima are not numerous enough. Actually, the thermodynamic transition 
is precisely due to the vanishing of the complexity at T s . Below T s the system is in 
the so-called 1-RSB phase (correspondingly the high temperature phase is called replica 
symmetric (RS)). For T s < T < T d , on the other hand, the complexity takes a finite 
nonzero value. In other words, between T s and T d , a.n exponential number of metastable 
states contributes to the free energy. In this temperature range, the system is already 
not ergodic since a configuration starting in one of these multiple states and evolving 
with say, Langevin dynamics remains trapped inside the initial state forever (i.e. on all 
timescales not diverging with the system size). 

When approaching T d from above, the majority of the stationary points of the TAP 
free energy, which are unstable above T d , become marginally stable at T d and stable 
below. One therefore expects a slowing down of the dynamics due to the rarefaction of 
descending directions in the free energy landscape [31] • In fact, the dynamics of some of 
these models, e.g. the p-spin spherical model, can be analyzed exactly. One can show, 
from the dynamical equation for the spin-spin time-dependent equilibrium correlation 
functions, that at T d an infinite plateau appears in the dynamical overlap q d (t) defined 
as: 




(1) 
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where <jj is the spin at site i. Remarkably, the integro-differential equations derived 
with the Mode Coupling Theory of Glasses for the density-density correlation functions 
reduce identically within the so-called schematic approximation [32] to the equation 
obeyed by q d (t). The solution of the equation for qd(t) (or the more complicated full 
MCT equations) [33] leads to a two-step relaxation for the correlation function: there 
is a first rapid decay from 1 towards a plateau value qEA, then a slow evolution around 
it (the f3 regime) and, eventually, a very slow decay from it (the a relaxation). The 
plateau value is called the Edwards-Anderson parameter, in analogy with spin glasses. 
For our present purpose, we will be interested in these two last regimes. 

At a given temperature T > T d , the dynamics in the /3-regime is described by 
power-laws. The approach to the plateau can be written as q d {t) ~ qEA + ct~ a , and the 
later departure from the plateau as qd(t) ~ q E A — c't b . The exponents a and b are model 
dependent, but satisfy the universal relation [31] 19]: 

r 2 (i + fr) = r 2 (i-q) 

r(l + 26) r(l-2a) ' 1 ' 

where T(x) = J °° t x ~ 1 e~ t dt is the Gamma function. In the a-regime, close to T d , 
qdif) verifies a scaling law, called time-temperature superposition in the structural glass 
literature: 

q d (t) ~ f(t/r a ) . (3) 

This scaling form allows one to superpose on a single curve the data for different values 
of time and temperature. A good fit of the function f(x) is obtained with a stretched 
exponential f(x) oc exp(— x 13 ). The a relaxation time r a diverges at the transition 
T = T d , as: 

r a oc (T - T d y , 

The dynamical correlation q d (t) is the order parameter of the dynamical MCT 
transition. Recently, especially in connection with the phenomenon of dynamical 
heterogeneity in glass-formers [35], there has been a lot of interest in the critical 
fluctuations and dynamical correlations associated with this transition. A central 
observable introduced in (36j |37] is the four-point susceptibility x^if), which measures 
the thermal fluctuations of the order parameter q d (t): 



XA® = N((q d (ty) - (q d {t)Y) (5) 

where, as usual, the brackets denote a thermal average and the over-line denotes an 
average over the quenched random couplings. For mean field glass models, or more 
generally for MCT (without conservation laws), one can show that the four-point 
correlation function becomes critical near the MCT transition temperature Td [38] [Tl] . 
In the /3-regime, and with e = (T — T d )/T d , 

X*(t) ^ —r /l(te^) > t~T P = (6) 
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and in the a-regime, 

X4 (t)~ ±f 2 (W) , t~r a , (7) 

where fi(x) and f2(x) are two scaling functions, with the following properties: fi(x) oc 
x a when i < 1 and scales like x b when x ^> 1; f'2{x) scales like x b for i < 1 and 
vanishes for large x (note that the large time limit of x^{t) is equal to the spin-glass 
susceptibility which is not critical at T d ). Using the definition of 7, and equation 
the numerical analysis of x±(f) ver Y close to Td allows, at least in principle, to extract 
all the MCT-exponents without having to know the value of qEA |38j . 

Finally, in [TUl [TT] it has been shown that MCT can be interpreted as a mean field 
approximation of a critical model. Following this point of view, one can compute, within 
this mean field theory, the spatial dynamical correlations and the associated diverging 
correlation length at the transition £(T). It can be established that £(T) oc 1/(T — Td) v 
with v = 1/4 (see also [39]). Furthermore, the upper critical dimension of the theory 
turns out to be d u = 6 for dynamics without exactly conserved variables [T31 UJ] . 

3. Numerical simulations of the Random Orthogonal Model: equilibration 
and static properties 

3.1. Definition of the model 

The ROM [18] is a fully-connected spin model with quenched disorder, defined by the 
Hamiltonian: 

H=~^a i J ij (T j , J= l OW , 

where A is a diagonal matrix whose elements are equal to ±1 and are drawn according 
to: 

p (\) = P 5(X - 1) + (1 - P )6(X + 1) , (8) 

and O is an orthogonal matrix distributed with respect to the Haar measure (that 
we generated numerically by using the NAG routine G05QAF), p is a real number 
p G [0, 1]. The will use later the notation u = {Jij} to denote a given instance of the 
disorder. The original model corresponds to p — 1/2. The normalisation is such that 
Jfj = 1/N 2 TtJ 2 = 1/N, Vp. A detailed analysis for arbitrary p can be found in [26J. 
Both the static and the dynamic temperatures depend on p (see [2S] for more details). 
In the following, p is set to 13/32 ~ 0.4. This gives us higher transition temperatures 
than for p — 1/2 with a good separation of T s and T^: T s ~ 0.102 and ~ 0.177 
respectively||] With this value of p, the annealed entropy vanishes at a temperature 
close to T s (For p > 1/2 the annealed entropy vanishes exactly at T s ). 

I There might be other another transition to a full RSB state at lower temperatures, but we will not 
be concerned by this possibility. 
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The static order parameter of the 1-RSB transition is the usual static overlap q 
between two replicas, i.e. two equilibrium configurations {cr} and {r} characterized by 
the same quenched disorder. Its probability distribution function can be written as: 

P(q)^(5(q-^Y1 (T ^- ( 9 ) 

i 

In the thermodynamic limit, P(q) = S(q) for T > T s , two different replicas have zero 
overlap with probability one. In the low temperature (1-RSB) phase T < T s , a second 
delta function peak centered at a value q\ > appears with a weight 1 — m. To wit, 
two different replicas may have a mutual overlap q = with probability m, and q\ with 
probability 1 — m. The Edwards-Anderson order parameter qEA(T) is equal to q\. Note 
that the shape of P{q) is only sensitive to the thermodynamics, and for T s < T < T^, 
P(q) = 5(q) like in the the RS phase. 

The ROM is strongly discontinuous: the value of q\ jumps sharply from to a 
value close to 1 at the static transition (see figure |6| below. This, together with the 
wide separation between T s and Td, are strong cases to use the ROM, as compared to the 
p-spin model [T7] or the Potts glass [20J |2H [22]- Both these models have unfortunately 
very close static and dynamic transition temperatures. 



3.2. Numerical method 

Let us start by giving some details about our simulations. We study systems with 
N = 32, 64, 128 and 256 spins. We thermalize the system using the parallel tempering 
optimized Monte Carlo procedure [HI HOI HH 112] , with a set of 100 temperatures in the 
range [0.08, 0.352] (AT = 0.002 for T < 0.2 and AT = 0.004 for T > 0.2), except for the 
largest system where a smaller set of 35 temperatures has been used, in order to save 
computer time (with 0.12 < T < 0.28 and AT = 0.008). These parameters have been 
chosen empirically and no claim is made that they are optimal. As usual, the program 
simulates the independent evolution of two clones/replicas, in order to compute the 
static overlap q. We perform 15 10 6 parallel tempering iterations (one iteration consists 
of one Metropolis sweep of all spins, followed by one tempering update cycle of all pairs 
of successive temperatures). The second half of the equilibration procedure is used 
to measure the static quantities. The results of this procedure are presented in the 



following subsections (see subsections 3.3 and 3.4). 

We then store the final hopefully well-equilibrated configurations. Due to the 
difficulty to equilibrate large systems, we restrict ourselves to T > T d to study the 
finite size dependence of the dynamics. We also restrict the number of temperatures 
to 25, equally distributed between 0.178 and 0.29 (AT = 0.008), for all system sizes. 
We then perform 5 10 5 pure Metropolis sweeps and measure qd(t), the overlap between 
the well equilibrated initial configuration, and the configuration at time t. We then 
perform 5 10 5 parallel tempering iterations in order to have a new, well de-correlated, 
starting point, and repeat the procedure 100 times. This gives us for each disorder 
samples N t h er = 200 thermally independent estimates of the dynamical overlap qd{t), 
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this is large enough to obtain a reliable estimate of the non-linear susceptibility \±. The 
number of disorder samples uj is = 500 for all sizes. 

Let us mention that we have also studied the model obtained by projecting the 
elements of the matrix J to ±1. The motivation was that in this case one could use an 
efficient multi-spin coding technique [33] . Unfortunately, the resulting model displays 
a very different physics, with oo-RSB breaking, and is accordingly not suitable for our 
purpose. This is an illustration of the fragility of the ROM model with respect to 
perturbations. 



3.3. Relaxation and equilibration tests 

Disordered systems are notoriously difficult to simulate, and it is crucial to ensure good 
quality sampling, and in particular good thermalization. There is unfortunately no 
full-proof heuristics for this purpose. 

The heuristics we use for this simulation is to check that the fluctuation-dissipation 
relation, relating the specific heat to the variance of the internal energy, is satisfied. 
This is a very stringent test (see e.g. [44]) and we have checked in the case of the SK 
model that it is consistent with other methods used in the literature. We define the 
ratio R(T) as: 

where e is the energy density. R(T) vanishes when the configurations are well sampled. 
Our data for R(T) (with e as measured in the second half of our thermalization runs) 
as a function of T can be found in figure [T] for system sizes N = 32 up to 256. The 
results are very satisfactory for all temperatures up to iV = 128. For the largest system 
size however, the sampling is clearly not good enough below T^, in spite of intensive 
numerical efforts. This is to be contrasted to the cases of the SK and Potts glass models 
where values of N up to a few thousands can be handled with the parallel tempering 
algorithm [351 120] . 

The modest efficiency of the parallel tempering algorithm when applied to the 
ROM can be directly observed by studying how the main thermodynamic observables 
reach their equilibrium values. Starting from a random initial configuration for the two 
clones (we take Oi = 1, Vi), we plot the instantaneous value of the considered observable, 
averaged over the disorder, in order to tame the fluctuations. The results for the internal 
energy (e(t)) and for the overlap between the two clones (q(t)) are given in figure [2J for 
T = 0.178 ~ T d . This figure suggest a power-law behaviour e{t) - oc t~ a3 in the 
large N limit, with a non uniform convergence. The smaller the value of N, the earlier 
in time do the data deviates from the power law behavior. For large t, the relaxation 
becomes strongly iV dependent, smaller systems relax much faster. 

A power-law relaxation of the energy at is in fact expected on theoretical 
ground. The value of the exponent was recently conjectured by A. Lefevre (private 
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Figure 1. Fluctuation-dissipation ratio R(T) (see equation 10 for the precise 
definition) as a function of the temperature, for N — 32 up to N = 256. We 
note that for all vales of N, R(T) fluctuates around zero, namely the systems are 
well equilibrated, except for the largest size (N = 256) where the system fall out of 
equilibrium below T c i ~ 0.177. 



communication) to be equal to a, the MCT exponent denned above. This is compatible 
with our finding since, as we shall find later, a ~ 0.35 for the ROM. 

The slow relaxation of the internal energy, even for T > Tj, illustrates why this is a 
hopeless task to simulate the ROM below (and near) Td for reasonable system sizes. The 
much faster relaxation of the overlap can be understood with the following argument: 
the two clones starts at the top of a very rugged energy landscape. After a few sweeps, 
they have started falling down in directions, or towards, traps whose probable overlap 
is zero. This leads to a fast decorrelation of q(t). Instead, in order to equilibrate the 
internal energy, one has to visit many different traps in order to have a good statistical 
sampling of all energy states, which is a much more slow process. 

Many reasons can be invoked to explain the poor performances of the parallel 
tempering algorithm applied to the ROM. The first one is that the temperature in 
the interesting region is in fact extremely low, leading to extremely small Metropolis 
and exchange acceptance rates. A way to gauge the smallness of the ROM transition 
temperature is to study the spin-glass susceptibility xsg — N(q 2 ) as a function of 
the temperature (see figure |3|. When T — > oo, xsg must tend to unity, but strong 
corrections are still present up to T « 10T d . This is expected, since the leading non 
trivial 1/T correction to xsg is the same in the ROM and SK models, because the 
variance of the is normalized to the same value in both cases. However, the critical 
temperature of the SK model is T^ K = 1, ten times larger than the critical temperature 
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Figure 2. Left: Relaxation of the disorder averaged internal energy from a random 
initial configuration, using the parallel tempering algorithm, for N = 64, 128 and 256 
dXT — 0.178, just above T^. The apparent power law behavior e(t) —e^ oc t~ 3 is only 
valid in the early time region. One notes that even (slightly) above T^, the relaxation is 
extremely slow and equilibration is quite difficult to achieve (one needs a couple of 10 5 
Monte Carlo sweeps in order to equilibrate the internal energy for N — 256). Right: 
Relaxation of the disorder averaged overlap between two clones from disordered initial 
configurations for N = 128 and 256 at T — 0.178. The relaxation is much faster (a 
few hundred Monte Carlo sweeps are enough to equilibrate the overlap). We used the 
same color code of figure 1. 



of the ROM. 

A deeper explanation relies on the physics of 1-RSB models, recalled in the above 
section. The existence of numerous metastable states in the range [T s ,Td] slows down 
dramatically the dynamics. Indeed, in order to have a decent sampling, one should 
explore a representative subset of all the metastable states. The complexity of the 
ROM has been computed in the thermodynamic limit in [16]. It shows a sharp jump 
from to a finite value below Tj, and thus the log of the number of states that must be 
explored jumps from to a number of order iV at It does not come as a surprise that 
the algorithm fails below even for moderate values of N. A precise understanding of 
why the ROM case is so much difficult than the Potts and p-spin cases is still lacking. 
A reasonable conjecture is that this is due to fact that the overlap value is so close to 
unity. Another system where this happens, and where the dynamics is indeed painfully 
slow, is the Bernasconi model (see [4T1 H8] ). 

We note, en passant, that applied to the model with binarized exchange couplings 
briefly mentioned above, the parallel tempering works brilliantly. 

3.4- Thermodynamics 

We show in figure [4] our data for the internal energy per spin, defined as u(T) = (H) /N 
as a function of the temperature, together with the theoretical result obtained in [26J. 
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Figure 3. Spin-glass susceptibility for the ROM as a function of the temperature, 
for N = 32 to 256. One observes (a) strong size-dependence effect below Td ~ 0.177 
and (b) an extremely slow convergence towards its asymptotic T — > oo value. In the 

(N—32) 

inset Xsg ^ s Pitted for high temperatures up to 35T S together with the asymptote 
X^sg ^ = 1- This illustrates how "cold" the ROM glassy phase is. We used the same 
color code of figure 1. 



For our special choice of p, one finds: 



u{T) = - a \T-\IT*-^T + A\ ,T>T S , (11) 

U (T) = - 0.47 , T <T S . (12) 

The numerical data for u(T) is qualitatively consistent with the infinite-volume 
analytical results, up to the finite size corrections, with the marked exception of the 
N = 256 data below T^, which are at odds with the rest of the picture. This is 
in agreement with our previous observation that the N = 256 systems are not at 
equilibrium at low temperature and should be discarded there. The temperature below 
which the fluctuation-dissipation relation is violated indeed roughly coincides with the 
one below which the numerical values of u(T) become manifestly wrong. 

After discarding the bad data, u{T) converges towards the predicted infinite- volume 
value, although with marked size-effects below the dynamical temperature. The finite 
size effects at T = 0.08 are roughly compatible with a 1/N behavior. A similar 
observation was made for the p-spin model in |45j . 

Following jlU] . we plot in figure ^ the coefficient A(T) defined as: 

A(T) = Ml^S , (13 ) 
(I 2 ) 
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Figure 4. Internal energy density u(T) as a function of T for N — 32 to TV = 128 
(top to bottom) and 256 (points that do not extend below T — 0.12), together with 
the analytical prediction in the infinite volume limit. 



that signals [50] the onset of the non self-averaging behavior of q 2 . In figure [5j we plot 
the usual Binder parameter, 

B(T) = ;(3-^| . (14) 
(? 2 > 




Figure 5. The parameter A (left) and the Binder parameter B(T) (right) as a function 
of T, iorN = 32 to 256. We used the same color code of the previous figures. 

For generic 1-RSB transitions, Picco et al. [H] have argued that these two 
coefficients are zero for T > T s (in the thermodynamic limit), non zero for T < T s 



Predictive power of MCT .-Numerical test and Finite size scaling for a mean field spin glass 14 

and in fact diverge at the static transition T — > T~. The fact that the Binder coefficient 
is negative when T — ► T~ is simply related to the appearance of a peak in P(q) for 
q 7^ 0. In the infinite volume limit, one has: 

P(q) = (l-m)8(q)+m5(q-qi), (15) 

with m < 1, and m — > when T — > T~ . Accordingly one has B(T) = 1/2(1 — 1/m) < 
in the 1-RSB region. 

Our data are in qualitative agreement with the above limiting behavior (see jlHJ [17] 
for a similar numerical analysis in the case of the p-spin model). Note that for a 1-RSB 
transition, the curves of B(T) as a function of T for various values of N do not cross at a 
universal point, at variance with usual phase transitions (including oo-RSB transitions). 

The overlap probability distribution P(q) is found to have the shape corresponding 
to a 1-RSB phase transition, with one peak around q = and, below T s , another peak 
around q — qi (qi is close to one in our specific case). With iV finite, both peaks have 
a non zero width, as always. It turns out that the q = peak is much broader that the 
peak at q = q\. Above T s , our data shows a spurious peak centered at q ~ 1. This peak 
however corresponds to an unstable thermodynamic phase and decays sharply with the 
size of the system (see figure [6]) as it should. As shown in figure (|6]), the peak centered 
at q ~ 1 becomes extremely sharp for T < T s . 
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Figure 6. Left: P(q) for T s <T = 0.154 < T d .. Right:P(g) for T = 0.078 < T s . Both 
for N = 32 to 128. We used the same color code of the previous figures. 



3.5. Thermodynamics: Conclusion 

We have thus been able to confirm numerically the main replica predictions for the 
ROM: the energy as a function of the temperature freezes at the static transition, and 
the order parameter is strongly discontinuous there. This could be related to the fact 
that the static transition temperature is very small, compared, for example, to the 
SK model. Correspondingly, it is very hard to equilibrate the system despite intensive 
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numerical efforts. Systems with N = 256 did not reach equilibrium below the dynamical 
transition. 

4. Numerical Simulations: dynamical behavior 

We now present our numerical study of the equilibrium dynamics of the ROM. On 
theoretical grounds, and in view of the above results on the statics of the model, the 
dynamics of the ROM should be described by the Mode-Coupling Theory, at least in a 
Landau sense, i.e. the MCT power laws are expected to be valid and the actual values of 
the exponent, although not universal, are constrained to verify equations [2] anc{4| see [9]) 
We compare our data with the predictions of MCT both for the two-point and four-point 
correlation functions. The results are very puzzling at first sight. We will show in the 
following sections that a detailed understanding of preasymptotic corrections and finite 
size effects is required in order to rationalize our numerical results. 

4-1. Dynamic scaling and comparison with MCT 

We first focus on the dynamical overlap qd(t) (defined in equation [I]). Our data show a 
plateau in qd(t), whose extension increases by lowering the temperature, see Fig [7} 




Figure 7. Left: Dynamical overlap qd(t) for the largest system (N = 256), and for 
various temperatures from T ~ T<j (upper curve) to T = 1.65 Td (The temperature 
step is such that AT/Td = 5%). Right: Test of Time-Temperature superposition for 
qd(t) for all sizes (N = 256 grey circle, 128 blue upper triangle, 64 red squares and 32 
black lower triangle) and various temperatures from T = 1.25 Td to T = 1.65 Td- One 
remarks clearly the decay of finite-size corrections and the convergence to a limiting 
curve. 



The value of the dynamical overlap on the plateau is close to the infinite-volume 
limit of the Edwards-Anderson parameter (qEA{Td) ~ 0.955), as it should. After a 
single Monte Carlo sweep, qd(t) has already decayed to the plateau value. Therefore we 
cannot observe the early /3-regime. This is an unfortunate drawback of our choice of the 
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ROM with parameter p = 13/32. We however do see the a regime in full glory. Our 
results are quite different from those obtained for the Potts model [20j [2TJ [22] , where 
no plateau was observed in qd(t) for systems with up to iV = 2560. 

In order to be quantitative, we define the time scale r a (T) (see equation [3]) as the 
time needed to reach the value qd(j a ) = 1/2. If Time- Temperature superposition (TTS) 
holds, the precise definition used is irrelevant. We have checked that the gd(t)'s plotted 
as a function of t/r a (T) approximately collapse on a unique scaling curve for the largest 
system sizes, but for temperatures not too close to Tj, see Fig [7} 

However, as T gets closer to T d , TTS appears to breakdown: instead of approaching 
a universal scaling curve, the qd{t)'s are more and more stretched as the system size 
increases. This is a priori surprising since TTS should work better close to the MCT 
transition. But we find that close to T^, finite size effects become important. This is 
also revealed by the behavior of r a (T) as a function of e = (T — Td)/Td, see Fig [8j A 
MCT power law fit, r a oc e -7 , accounts reasonably well for the regime e > 0.2 where 
finite size effects are small, and yields 7 ~ 2.1. MCT also makes detailed prediction 
about the form of the scaled relaxation function, as we discussed in the introductory 
sections. Unfortunately the power-law predicted in the early j3 regime is inaccessible 
because the plateau is too close to unity. However, we can check the (von Schweidler) 
power-law associated to the late (3 regime. In order to do so in a way that does not 
require a precise determination of the plateau value qEA, we plot in figure [8] — dqd/dhit 
for our largest system size N = 256 at the lowest temperature at which we do not have 
substantial finite size effects. According to MCT, this quantity should increase as t b 
when < ( < r a . From the power law fit shown in figure [8] we determine b ~ 0.62. 
But within MCT the values of b and 7 are not independent (see equations (2]4) above). 



The value of b corresponding to 7 ~ 2.1 is found to be &mct ~ 0.75 which is distinctly 
too large compared to our data (see Fig [8]). This is a second puzzling result since MCT 
predictions are expected to apply to the ROM dynamics close to T d . 

Other puzzling features emerge from the analysis of the dynamic susceptibility 
% 4 (t), as defined by equation A plot of the peak value xl as a function of e for our 
largest system size shows (in the regime without finite size effects) a power law behavior 
compatible with the MCT prediction xl oc 1/e (see figure [9]). However, figure [9] shows 
that there is no collapse of X4(0/xl plotted as a function of t/r 4 , where T4 is such 
that xl = Xi{t — r 4|§} a t variance with the MCT prediction for t < r 4 Before the 
peak, Xiif) appears to grow as t bi with 64 « 0.85. Such a power-law increase is again 
predicted by MCT, but one should find x±{P) °c —dqd/dlnt [101 HI], i.e. 64 = b w 0.62. 



In fact, MCT also predicts that x^if) oc —dq d /dT [14]. We show in figure 10 these three 
different quantities in log-log plot. Although —dqd(t)/dT and —dqd/d\iat are indeed 
similar, x^if) does not conform to expectations. 

In conclusion: although some of the MCT predictions are quantitatively obeyed, 

§ We checked that T4 and r a are approximatively proportional. 

|| For large times scaling is not expected, since in this limit Xi(t) ~> Xsg an d this has a finite nonzero 
limit as N — > 00 and T — > Td, contrary to xt that is expected to diverge. 
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Figure 8. Left: Relaxation times r a versus e = (T — Tdj/Td, for different system sizes 
(same color code of previous figures). For e > 0.2, i.e. T > 1.2 Td, an TV-independent 
regime is observed, with r a oc e~ 7 and 7 ~ 2.1, over two time decades. Right: 
Determination of the 6-exponent of the von Schwcidler law through the derivative 
of the overlap, -dq d /d\nt oc t b . (Here N = 256, T = 0.226). The two straight lines 
correspond to b = 0.62 and buCT = 0.75 respectively. 




Figure 9. Left: Maximum of \4 versus e. The MCT prediction would be xt ^ V e - 
Right: Scaling plot for Xi (here N=256). The straight line corresponds to a power law 
with an exponent 0.9. 



other important ones are clearly violated. The solution of this conundrum is that the 
values of e used above are actually not small enough to be in the asymptotic regime 
where MCT predictions hold. As we shall show in the next section, these predictions 
are only valid in a surprisingly small region close to the transition. So why not work 
closer to T^? The next problem we will have to deal with (section [5]) is finite size effects, 
that become large close to T^. Therefore, only after a very careful finite size analysis 
can one conclude on the compatibility between the MCT predictions and the numerical 
behavior of a model that is in principle exactly described by MCT! We will show how 
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difficult this program turns out to be for the ROM. This sheds considerable doubt on 
the precise, quantitative comparison between experimental data and MCT, since these 
problems should show up in these well. 




0.01 1 1 

0.01 0.1 1 10 

t/T 



Figure 10. Comparison of three susceptibilities: \4 (circles), —dq^/dlnt (squares) 
and —dqd/dT (triangles). Here, N = 256 and T = 0.226. Each susceptibility has been 
rescaled with its own \* and t* . All three quantities should grow with the same initial 
power-law behavior within MCT. This is clearly not the case for X4(t). 



4-2. MCT critical properties and preasymptotic corrections 

At this stage, it is important to have a reliable reference point to which we can compare 
our numerical results. For this we choose to study in detail the Leutheusser integro- 
differential equation for the correlation function jJT] , which comes out of the schematic 
version of MCT with a so-called quadratic kernel: 

$(t) +m(t) +$(*) +4A f rfr$ 2 (r)$(t - r) = 0, $(0) = 1 (16) 

Jo 

In the above equation, $(£) is the correlation function and plays the role of qaif) above, 
and A is the coupling constant that measures the strength of the feedback effects at the 
heart of the MCT transition. Remarkably, the equation above is also the one governing 
the evolution of the correlation function for the p = 3 mean field disordered p-spin model 
[28] . It can be analyzed mathematically, and all the results quoted in section [2] can be 
shown to hold exactly in the limit A — > A^ = 1. In particular, the model is ergodic for 
A < 1, where lim^oo $(£) = 0, and develops power-law regimes with exponents a and b 
given by (equation [2]): 

r 2 (i + &) _ r 2 (i-q) _ i 

r(l + 26) T(l-2a) 2' 1 ' 
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leading to a « 0.315, b — 1, 7 « 1.765. 

The Leutheusser equation can be solved numerically for arbitrary large values of t, 
using for example the algorithm of [52]. In what follows we fix Q = 1, and neglect the 
$(£) term, as usually done. In order to compare directly with the ROM data above, 
we have computed $(t) for values of A that are at the same relative distances from the 
critical point as our ROM data, for temperatures T = 0.178, 0.186, . . ., and 0.250. (We 
recall that T d = 0.177 for the ROM.) 




le-06 0.0001 0.01 1 100 0.1 1 10 



t/T a t/T a 



Figure 11. Leutheusser model. Left: Scaling plot of Q(t) as a function of t/r a , 
at the same relative distances from the MCT transition as for our ROM data above 
T d , namely A = T d /T, for T = 0.178, 0.186, . . . , 0.250. The a time scale is given by 
T a = (1 - A) -1765 . Right: Scaling plot of -d$(t)/dlnt. The expected late (3 behavior 
t b with b = 1 is shown for comparison. The straight line represents y = t/(5r a ). 



Figure 11 shows a TTS plot of using the exact value of 7, i.e. using r a = e~ 1-765 , 
where e = (1 — A) = (T — T^/T^. This figure shows apparent scaling in the late (3 (von 
Schweidler) regime, with only small scaling violation. However, only the relaxation 
corresponding to the value of A closest to unity (e ~ 0.05) reveals the expected two-step 
relaxation with a nontrivial plateau region! 



More revealing is a log-log plot of — /dint, shown in figure 11 together with the 
expected theoretical behavior — d$>(t)/d\nt oc t b with 6=1 (The derivative is computed 
by plain finite difference). We now see very strong scaling violations before the peak, 
and an apparent value of b that is significantly below 1, even for the curve closest to the 
critical point. The conclusion is that while the value of 7 extracted from a TTS plot 
of $(t) is reasonable, the value of b that one can extract from $(t) 5% away from the 
critical point is grossly underestimated. This is similar to our observations above for 
the ROM. 

Let us now turn to the non linear susceptibilities x±{t) an d Xr(t) = d&(t)/d\. 
Within MCT, both quantities have the same scaling behavior. In particular one expects 
that Xi(t) = t~ l F{t/T a ), with F{ x) oc x b , in the late (3 regime. Numerically, Xrit) is 
easy to obtain from the value of $(£) for different values of A. The case of x&{t) ls ^ ess 
straightforward. It turns out [12] that a dynamical susceptibility, proxy of Xzi(t), can 



Predictive power of MCT .-Numerical test and Finite size scaling for a mean field spin glass 20 



be computed for the spherical p-spin model, which as recalled before is characterized 
by a dynamical equation for the correlation function identical to the the one of the 
Leutheusser model. A perturbed-time dependent Hamiltonian H v is introduced: 

H V = H- md {t) (18) 
where H is the usual p-spin Hamiltonian, and q<i{t) is the overlap between the spin 
configurations at times t and t = 0. Then Xiif) = d(qd(t)) / di] . For a given value of rj, 

rj=0 

one is led to a set of two integro differential equations that can be solved numerically [53J . 
The estimate of X4,{t) follows from a careful extrapolation to rj = 0. Figures 12 shows 




Figure 12. Left: Scaling plot of xr(t) = d$(t)/d\. The values of A are the 
same as in figures |11| with the same color code. The expected scaling behavior 
d$(t)/dln e cx (t/8T a ) b with b = 1 is shown for comparison. Right: Scaling plot of 
X4(t) for the spherical 3-spin model, with again the same values of e and the same 
color code. The straight line corresponds to the theoretical prediction \4 oc e~ 1 (te y ) b . 

scaling plots for xrif) and x^if)- Although an approximate scaling is observed close to 
the peak, only the Xt(^) curve closest to the transition gives a hint of the correct value 
of the exponent b. The scaling violations are non monotonous and can fool the reader 
into seeing scaling with some b less than the correct value. The estimate for b from 
Xi(t) is systematically larger, and closer to the true value. This again is similar to our 
numerical observations for the ROM. 

In order to observe the asymptotic MCT scaling predictions, one must work much 
closer to the transition. For example, the expected linear regime of — d§(t)/dlnt oc 
(t/r a ) only appears very slowly as A — > 1, and is well developed only for e < 10~ 4 . 



Figure 13 shows a similar behavior for d&(t)/d\. The linear region appears only when 



e < 10 . Finally figure 13 shows the even slower approach to scaling of x^if) (beware 



however that e is here limited to 10 6 ) 
4-3. Dynamics: Conclusion 



The above results show that the true asymptotic regime of MCT is unusually narrow. 
This should be remembered when comparing numerical (or experimental) data with 
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Figure 13. Same as figure [12J but going deeply inside the scaling region, namely for 
e = 1CT 1 , 1CT 2 , . . . , 1CT 8 for the left wing figure, e = 1CT 1 , 1CT 2 , . . . 1CT 6 for the right 



wing figure. The straight line arc the same as in figure 12 



MCT predictions. These data are usually plagued with noise and possible finite size 
effects, on top of the strong scaling violations that appear even in the best case situation 
studied in this section. On the other hand, some useful conclusions emerge, which allows 
us to make sense of our data on the ROM dynamics in the preasymptotic regime: (i) 
an approximate TTS holds for the late j3 regime, with the correct value of exponent 
7; (ii) the exponent b extracted from the time dependence of the correlation function 
underestimates the true value; (iii) for the dynamical susceptibility the scaling is 
acceptable around the peak with the predicted divergence e _1 , with a value of b in the 
correct range. From these considerations, we conclude that the correct values of b for 
the ROM should be around b « 0.8, and the corresponding value of 7 close to 2, and 
a »s 0.35. A confirmation of these values should come from studying the dynamics closer 
to Trf. However, finite size corrections become important there and we now turn to the 
study of these effects. 

5. Finite-Size Scaling: More surprises 

In the previous section, we have shown that the critical behavior of the ROM dynamics 
in the regime where finite size corrections are small is polluted by strong preasymptotic 
effects. In order to get rid of those one should simulate very large systems very close 
to Td, which is alas not possible since the equilibration time also becomes very large. 
The hope would be to use finite size scaling (FSS) to extract the interesting asymptotic 
behavior. 

5.1. Naive theory and comparison with numerical data 

The MCT predictions are modified for finite but large system sizes. One expects in 
particular that activated effects, absent for infinitely large systems, start playing a role 
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for finite systems close to T^. 

As we have recalled above, it was recently recognized that MCT is a mean field 
(Landau) theory characterized by a diverging length scale £ oc e" 1 / 4 [TOl ITT] and the 
corresponding upper critical dimension is d u = 6. Assuming that the field theoretical 
analysis of [Ml ES] applies also to this dynamical transition ffij] we expect that finite size 
scaling holds for MCT above the upper critical dimension d u where the proper scaling 
variable is not L/£ = N l / d e u but rather N 1 l du e v [SB]- The fully connected ROM is 
obviously above the upper critical dimension and the relaxation time r a for a finite 
system should therefore take the following scaling form: 

r a (T, N) = N^kF(N^e), vd u = 3/2 (19) 

When N — > oo, all N dependence should disappear and the MCT divergence of 
equation^ must be recovered. This means that the scaling function T must behave as 
T[x) oc x -7 when x ^> 1. 



We analyze our numerical results on r a (T, N) using the above FSS form, see Fig 14 



A good collapse of the different curves can indeed be obtained using equation ( 19 ) with 
7 = 2.1, but we need to use the value d u u = 2 instead of the expected value 3/2. 




Figure 14. Left: Scaling plot of r a (T,N) using equation (19 1 with 7 = 2.1 and 
d u v — 2. Right: Scaling plot of the maximum of x 4 , using equation (20 1 with d u v = 2 
and s = 1. 



If we now turn to the four-point susceptibility X4, another confusing result is 
obtained. While the peak location has the same finite size scaling expected, 
the FSS of the peak height should read: 

XI(T,N) = N s G(N^e). (20) 

Using the above effective value vd u = 2, the exponent s should be such that for large 
N, Xa diverges as e _1 independently of N. This fixes s = 1/2, at variance with our 

Note that although a direct field theoretical analysis of FSS for the MCT dynamical transition seems 
very difficult, the scaling MCT exponents are related to the ones obtained from the replica theory. 
Naively, the analysis of [S31 [SS] is expected to hold for the replica field theory. 
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numerical data that suggests xt <x N s with s ~ 1, see figure 14 Note that despite the 
uncertainties on the exact value of s, we clearly find that s > 1/2 (note also the decay 
of the scaling function Q{x) at large arguments). 

We therefore find that finite size scaling appears to work for the ROM transition 
but not in the way expected, at least naively. In particular, the scaling variable \/N e is 
at odds with the value of the upper critical dimension for the ROM dynamics (without 
explicitly conserved variables) and the value of the exponent v. 

The next sections are dedicated to a detailed discussion and explanation of these 
puzzling results. As we shall show there is no contradiction whatsoever and the origin 
of this strange FSS are sample to sample fluctuations of the critical temperature T d . 

5.2. Harris criterion and random critical points 

It is well known that the stability of pure critical points with respect to weak disorder pTTj 
is governed by the Harris criterion: near a second order phase transition in dimension 
d, the bond disorder is irrelevant if the specific heat exponent a pure = 2 — du pure is 
negative, where z/ p?tre is the correlation length exponent of the pure system. Then the 
critical exponents of the disordered system are the same as the ones of the pure system. 

If a pure > 0, disorder becomes relevant and the system is driven towards a so called 
random fixed point characterized by a new correlation length exponent vr satisfying the 
general bound v R > 2/d. In the last twenty years, important progress [581 EH EQ] has 
been made in the understanding of finite size properties of random critical points. For 
our purpose we only recall that to each realization of disorder (u), one can associate 
a pseudo-critical temperature Td(u,N), defined for instance as the temperature where 
the relevant susceptibility is maximum. The disorder averaged pseudo-critical critical 
temperature = T d (u>, N) converges towards its infinite size limit as: — oc 
N~ l l dUR . The width AT^(A r ) of the distribution of the pseudo-critical temperatures 
Td(u), N) then depends on the nature of the critical point. If the disorder is irrelevant, 
ATd(N) scales trivially like iV -1 / 2 , but like N~ l / dVR if the disorder is relevant. 

All previous studies of FSS for disordered systems mentioned above have focused 
on models below their upper critical dimension (see for example |61J). Instead the ROM 
is clearly above its upper critical dimension. As a consequence, it is not obvious to 



deduce its FSS behavior from previous works. In Appendix A we discuss in detail the 
subtleties of the Harris criterion above the upper critical dimension. 

From a phenomenological point of view, one expects that general features of random 
fixed points should still occur: in particular one can define a sample-dependent pseudo- 
critical dynamical temperature Ta(uj,N) = Td + 5T(u,N). A hand-waving argument 
to understand the origin of these fluctuations is to consider the TAP equations for the 
ROM. The high temperature expansion leading to the TAP equations for the ROM has 
disorder-dependent corrections of order N~ l l 2 . These have a dramatic effect on FSS 
since these corrections are much larger than the expected FSS thermal window N~ 2 / 3 of 
the MCT-transition. We will therefore assume, and justify later on, that for each sample 
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T d (C4N) T d (£N) 

Figure 15. Sketch of the sample to sample fluctuations of a susceptibility in situations 
where the critical temperature fluctuates on scales N~ x / 2 , much larger than the FSS 
thermal window 7V~ 2 / 3 . Here lo and £ are two disorder configurations. 



the FSS window is indeed of the order iV~ 2 / 3 around a random critical temperature 
that has disorder fluctuations of the order N^ 1 ^ 2 . As a consequence, FSS for disordered 
averaged observables are dominated by the fluctuations of the critical temperature that 
wash out the much sharper TV -2 / 3 FSS thermal window [j (see figure 15 for a cartoon 



representation). Note that a similar situation for FSS of random first order transition 
has been discussed by D.S. Fisher in 



5.3. Random critical temperatures and modified FSS 

To make the above statements more precise, let us consider the generic example of 
some thermodynamic observable O. One would like to compute O. We assume that 
averaging over the disorder is equivalent to averaging over the distribution of critical 
temperatures, p(T d (u>, N)). For the sake of simplicity, p(T d (uj, N)) is taken to be a well 
behaved distribution with width of order iV _1//2 centered on Tj, the true asymptotic 
dynamical temperature Moreover, we posit that for each sample, some kind of FSS 
holds, in the sense that: 

0(T, N) = iV 2C/3 J^(iV 2/3 [T - T d (u, N))) (21) 

where ( is a certain exponent that depends on the particular observable, and the 
scaling function T is, at least to leading order, sample-independent. We assume that 
T{x) is regular for small arguments, and cx ^4± | | for x — > ±oo, such that 

0(T, N — > oo) oc |e|~^ independently of N. It is possible to justify all these assumptions 
within a simple toy model, the weakly disordered version of the Blume-Capel model, 

+ This cannot happen below the upper critical dimension due to the Harris criterion. 

* There will also be finite size corrections to the center of the distribution, but these are expected to 

be subleading to N~ x / 2 . 
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that displays exactly the unusual FSS discussed in this section. We warmly invite the 
reader to examine |Appendix B for more details. 



Now, writing T d (u, N) = Td — yN x l 2 , disorder averaging is obtained by computing: 

Q(T) = N 2<:/3 j dy p(y)F(N 2/3 e + N 1/6 y) (22) 

where e = T — T d . The analysis of the above integral in the large iV limit requires to 
distinguish two cases: £ < 1 and £ > 1 (with further logarithmic terms when ( = 1). 



• When ( < 1, one can set e = u/y N and take N 3> 1 in the above integral, to get 
to leading order: 

0(T) = N C/2 F(VN[T -T d ]), (23) 

A 00 , A+p(y) r u , An(y) , . 

T{u) = / dy + / dy 24 

J-u {u + y^ J-oo \u + y\^ 

Since ( < 1, the integral defining T is always convergent when u — > 0, leading to 

a well behaved scaling function. Therefore, in this case the usual FSS strategy is 

valid, with a scaling variable dominated by the fluctuations of critical temperature: 
AT 2 /3 _> jyi/2_ 

• When £ > 1, on the other hand, the relevant change of variable is y = z/N 1 ^ 6 — u. 
Again to leading order, this gives: 

0{T) w AN Vj, (-v/iV[T - T d ]) (25) 

where A = J dzT{z) is a convergent integral thanks to the rapid decay of T{z) for 
large z, leading to a finite multiplicative constant. In this case, FSS is drastically 
altered by the sample to sample fluctuations of the dynamical temperature; the 
decay of the scaling function is related to the one of the distribution of critical 
temperatures, and the exponent (4£ — l)/6 is unusual. 

The above analysis can be extended to the case of "assymetric observables" , where 
the power-law decay of Fix) is different when x — > +oo and x — > — oo. Most of the 



interesting ROM observables turn out to be of that type around T d , see section 5.5, 5.6 



5.4- FSS for single samples 

Checking these assumptions numerically is tricky for the ROM. If there was a quantity 
(susceptibility-like) with a sharp peak around T d then it would be simple: one would 
just rescale for each sample this quantity around its peak and verify whether the usual 



finite size scaling holds, as suggested by figure [L5j Unfortunately, no such quantity 
exists for the ROM. As we have seen, xt(T) is a monotonously decreasing function of 
temperature. In order to check the usual FSS, one should shift horizontally xt{T) for 
each samples around its own effective dynamical transition temperatures T d (u,N). In 
order to determine this effective critical temperature we focus on the sample to sample 



fluctuations of the relaxation time (see left panel of figure 16). We assume that the 



relaxation time is uniquely determined — at fixed — by the distance T — T d (u,N), 
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Figure 16. Left: r a is plotted as a function of the temperature for 5 different disorder 
samples (N = 32 with Nther = 120 only). There are clear strong sample to sample 
fluctuations of the relaxation time. Close to Td, the five curves can be superposed by 
a suitable horizontal shift, namely one has the relation t(uj,N) = C(T — Td(u>,N)). 
Right: Scaling plot of the relaxation time r a using rescaled sample dependent dynamic 
temperature. Here arc plotted the fastest, the slowest and a typical sample for all sizes 
(N=256 corresponds to grey circle, N=128 to upper blue triangle, N=64 to red squares 
and N=32 to lower black triangle). We recovers the standard FSS exponents by using 
Td{oJ,N) instead of the thermodynamic Td 



i.e. t(u,T) = C(T — T d (u,N)), where C is a certain function. Thus, choosing 
a certain reference relaxation time r* for a given sample, one fixes the difference 
T*(u,N) -T d (u,N) to a sample independent value 6* (with t(lu, T*(lu, N)) = r*) . 
Therefore, by averaging, 



5* = T*(u, N) - T d (u, N) = T*(lu, N) - T d (u, N) = T*(lu, N)-T d + 0(N 



-2/3A 



where we have assumed that the N dependent correction to the average critical 
temperature is 0(iV~ 2 / 3 ), that only introduces a shift in the final scaling variable. 
The above equation allows one to determine the sample dependent shift of critical 
temperature, T d {uj,N) — T d , as T*(u,N) — T*(u,N). This procedure does not require 
the knowledge of the functional form of the relaxation time. Once this shift is know, 
one can rescale the temperature axis in a sample dependent way, and test FSS sample 
by sample. 



The results are given in figures 16 and 17 Our statistics is quite limited 



because of the CPU-time consumption of such simulations. The modest number of 
thermal configuration does not allow us to extract the whole probability distribution 
for T d (u,N). Despite these limitations, the results are in perfect agreement with 
our expectations. After the above temperature rescaling, the best collapse of the 
relaxation time data is obtained with the naive finite-size scaling (N 2 / 3 (T — T d (cu,N))) 
and not y/N(T — T d (u,N)). A similar behavior is obtained for xl- The dynamical 
susceptibility divergence is now compatible with an e _1 behavior derived analytically 
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Figure 17. Scaling plot of the maximum of Xi with sample dependent transition 
temperatures Td (same convention as figure |l6| ). The behavior of xtiT) is now 
compatible with the MCT-predictions, x\{T) oc e _1 and with the standard finite-size 
scaling exponents. 



for MCT transition. 

5. 5. Anomalous FSS for the dynamical susceptibility 

Now we are confident in the validity of our analysis in terms of sample dependent 
transition temperatures, we come back on the anomalous FSS properties of the 



dynamical susceptibility that we observed numerically in section |5.1| above. We must 
first guess the shape of the sample dependent FSS for the peak susceptibility. When 
T > Td(u,N), we expect a divergence of xt(T) as e -1 , as discussed in the section 
above. So the exponent ( + corresponding to this region is C+ = 1- O n the other 
hand, below the transition one expects that the variance of the dynamical overlap will 
be of order N, due to activated dynamics that makes the system hop between states 
with zero mutual overlap (see the above remark on the fast relaxation of qd{t) in the 
ROM). Matching the requirement xl{T) oc N with the finite-size scaling form valid near 
T d (u,N), xt(T) ~ N 2 ^^(N 2 ^(T - T d (u,N))) leads to F(x) ~ for x -> -oo, or 



C_ = —1/2. Extending the analysis of section 5.3 to this strongly asymmetric case where 
C + 7^ we find that the average behavior is dominated by the left tail of T(x), finally 
leading to: 

xl(T)~N s g{VN{T-T d )), (26) 
where the anomalous exponent s is equal to 3/4 and Q is a certain scaling function. We 



therefore qualitatively understand the anomalous FSS result obtained in section 5.1, in 



Predictive power of MCT .-Numerical test and Finite size scaling for a mean field spin glass 28 



particular the fact that s is larger than the naive value s = 1/2. 
5.6. Relaxation time: some conjectures 

We finally turn to the relaxation time r a , which, as we already know, is very strongly 
sample dependent. The guess for the sample dependent FSS of the relaxation time must 
now account for the e~ 7 divergence for T > Td(u,N), and the activated dynamics for 
T < Td(u}, N). If we assume that the relevant energy barrier scales with the system size 
as N^, one is led to the following ansatz: 

t(uj,N) = N^F{{T -T d (w,N))N*), with (27) 
F{x) oc exp(— C(oj)\x\~) x — > — oo (28) 
J-{x) oc x~ 7 x — > oo, (29) 

where C(u) is a sample dependent constant that accounts for a possible sample 
dependence of the scaled barrier height. In the case of the ROM, it is reasonable 
to expect that if) = 1 [531 EU ESI [TT| . The above equation suggests that any observable 
governed by low enough moments of r(ou, N) will correspond to negative values of the 



exponent £ in the analysis of section 5.3 above. Since our definition of the average 
relaxation time r a is such that qd(t = r a ,T) = 1/2, this quantity is dominated by typical 
samples and we expect standard FSS with scaling variable \f~N~e, as indeed found in in 
section 15.11 

We however expect very different results for quantities sensitive to large relaxation 
times, dominated by rare samples. For example, the long time asymptotics of qd(t,T) 
is dominated by particularly "cold" samples. Neglecting the fluctuations of C(u), and 
assuming a Gaussian distribution of critical temperatures, we find: 



dy 



Qd(t,T) ~i^oo / — - exp 



\plixo 



2 



3± 

-IX ••■ ( 



?7 - C N^' 1 / 2 (y-u)^ 



(30) 



with u = yNe. Evaluating this integral by steepest descent, one finds that the 
asymptotic relaxation regime is, to leading logarithmic order: 



In q d (t, T) oc t _ 



hit 



4 



3i> 1 

(V> > (31) 



This decay is far slower than a stretched exponential, and gives a rationale to explain the 
observed slowing down of the late relaxation of qa, and is consistent with the behavior 
seen in the figures [7j In particular, we expect this slowing down due to cold samples to 
become dominant close to T^. 

6. Conclusions 



The aim of this study was to test numerically the predictions of MCT in a best case 
situation, namely for a model with an exact MCT transition, and analyze its finite size 
scaling behavior. 
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We chose the fully connected Random Orthogonal model, with a choice of parameter 
such that the dynamic (MCT) transition temperature is well separated from the 
static (Kauzmann) transition. We have first compared the theoretical predictions for the 
static (thermodynamic) properties of the model with our numerical results. Although 
we are not able to equilibrate large systems below T^, we find a good overall agreement. 
The transition temperature for the ROM is very low compared to the scale of the 
interactions; this implies that the transition is very strongly discontinuous, with an 
Edwards- Anderson order parameter very close to unity as soon as T < Tj. 

We have then studied the equilibrium dynamics of the model, focusing on the time 
correlation function and the four-point dynamical susceptibility, which measures the 
strength of dynamical heterogeneities. When comparing our numerical results to the 
predictions of MCT, we find that while some of these predictions are quantitatively 
obeyed (like approximate Time Temperature Superposition in the late (3 regime), other 
important ones are clearly violated, with inconsistent values of the MCT exponents. 
This is due to strong pre-asymptotic effects. Indeed, we have shown that the asymptotic 
MCT predictions are only valid inside an unusually narrow sliver around Td, thereby 
explaining these quantitative discrepancies and allowing one to get rough estimates of 
the MCT exponents for the ROM: b ~ 0.8 and 7^2. Working closer to Td to get 
rid of these strong preasymptotic corrections is hampered by equally strong finite size 
corrections. On that front, more surprises emerge: we find that the usual Finite Size 
Scaling (FSS) fails to account for our data, a result that we rationalize in terms of 
strong sample to sample fluctuations of the critical temperature. We have developed 
a phenomenological theory for FSS in the presence of these strong fluctuations. This 
modified form of FSS accounts well for our results; we also show that naive FSS works for 
individual samples. En passant, we have also developed new arguments to understand 
FSS in disordered systems above their upper critical dimension (see Appendices). 

The compatibility between the MCT predictions and the numerical behavior of 
a model that is in principle exactly described by MCT turned out to be extremely 
difficult to establish quantitatively, partly because of the impossibility to equilibrate 
large systems. The situation is expected to be worse when dealing with experimental 
data for which the critical temperature is blurred by non-mean field effects. In this case, 
quantitative comparison with MCT requires extreme care, to say the least. Our results 
show that some predictions appear to be more robust than others and this provides 
some guidance when dealing with application of MCT to experimental or numerical 
data. Indeed, we notice that the kind of violations of MCT predictions found in our 
study resemble very much what found in real liquids, see e.g. the difference between the 
time evolution of %4 and xt in [66J. 

On a different front, our results maybe relevant for FSS studies of super-cooled 
liquids [67J US]- In this case it has been shown that the usual theory valid for 
pure systems fails in account the finite size scaling behavior [15J. Although several 
justification can be put forward, in particular that the correlation length is not much 
larger than the microscopic length, our results suggest that new phenomena might 
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be at play. In fact, if the dynamically self-induced disorder present in super-cooled 
liquids plays somehow the role of the quenched disorder present for the ROM, as often 
proposed, then strong disorder fluctuations could lead to violation of usual FSS. Results 
qualitatively similar to the one reported in [15] (such as e.g. figure[3]) are indeed expected 
for the ROM. It would be certainly worth to pursuing further the comparative study of 
FSS in the ROM and real liquids. 

Finally, the study of a finite-range ROM that replicates the phenomenology of 
finite dimensional supercooled liquids is an interesting project that we are currently 
pursuing, in particular to test the predictions of the Random First Order Theory on a 
physical model that is as close as possible to its theoretical idealization. Other directions 
worth investigating numerically include a better understanding of the low temperature 
activated dynamics, which is probably only accessible in the aging regime. 
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Appendix A. FSS and Harris criterion for disordered systems above their 
upper critical dimension 

In the following we shall consider the role of disorder on finite size scaling above the 
upper critical dimension. This will lead us to formulate two different Harris criteria: one 
valid for the FSS and the other for the critical exponents, obtained from susceptibilities 
in the thermodynamic limit. 

Let us consider a pure system that is perturbed by the addition of small quenched 
disorder. As usual, we will focus on a disorder that couples to the energy, e.g. random 
couplings. As a consequence, samples of size N will be characterized by fluctuations of 
the critical temperature 8T C (N). 

In order to understand whether disorder affects the FSS behavior one has to 
compare the above fluctuations with the FSS window. Above the the upper critical 
dimension FSS is subtle [68]: the scaling variable is eN 1 ^ AuVu > where v u is the mean 
field exponent, e = (T — T c )/T c and d u is the upper critical dimension, see jUH] for a 
numerical check in five dimension for the Ising model. This means that properties of a 
pure finite systems depart from the ones expected in the thermodynamic limit when the 
distance from the critical temperature becomes smaller than j /V _1// ^ lll/ "^. Assuming that 
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on scales such that eN' 1 ^^ oc 0(1) the fluctuations of the critical temperature are of 
the order 1/y/N, one finds that for d u v u < 2 and close enough to the critical point, the 
addition of a small quenched disorder will therefore make fluctuates T c on a scale much 
larger than the FSS window of the pure system implying that the FSS behavior will be 
drastically affected by adding an infinitesimal disorder. 

This provides a generalization of the Harris criterion for FSS properties of systems 
above their upper critical dimension. As we shall show, contrary to what happens 
below d u , a different Harris criterion establishes when the disorder changes the critical 
properties of an infinite system. 

In order to investigate whether an infinitesimal disorder affects the critical 
properties, let's focus on the critical behavior of a generic local observable O x , e.g. 
the average local energy in x. A simple way to establish the Harris criterion consists 
in studying the perturbation induced by the disorder to the critical behavior. If the 
correction, no matter how small it is, ends to be the dominant contribution close to Tq 
this means that the disorder is relevant. Calling J y the random coupling at site y one 
obtains that the corrections due to the disorder are: 

y " 

This is a random variable whose typical value is given by 

2 

V (A.l) 



501 



X 



\ 



y 



'y 



where rj is the very small variance of the random couplings, J x J y = r\ 5 x ,y 

Below the upper critical dimension, just by scaling or using more refined techniques 
[68J, one knows that for the pure critical systems dO x /dJ y oc l/\x — y\ d ~ a l u where a — 1 
is the exponent characterizing the singular part of O x , which is \T — T c |~ a+1 . As a 
consequence one finds that the disorder fluctuations scale as ^^- d / 2 + a / u — <q€ dv / 2 ~ a . 
These will become dominant with respect to the pure critical behavior \T — T c \~ a+1 
when dv/2 < 1, no matter how small is 77, if one is close enough to the critical point. 

This is the standard Harris criterion. What does it change above the upper critical 
dimension? Actually, the previous derivation can be repeated identically. The only step 
where we used that d < d u is the assumption on the power law behavior of the response 
function dO x /d.J y . Above d u , one could just use the mean field critical power law 
behavior. This already would suggest that the Harris criterion for the critical properties 
is different from the one for FSS. However, the analysis is tricky because above d u one 
finds that sub leading corrections to the critical behavior dominate the sum in ( A.l ). Let 
us consider for instance the 4 field theory describing the Ising ferromagnetic transition 
and let us take O x = (<fi x ). We consider that the random couplings lead to a fluctuating 
mass in the field theory, i.e. the disorder couples directly to O x . In this case the response 
function dO x /dJ x+r reads below T c and for 1 <C r <C £: 
c d 

+ (t-t; 



r 2d-4 v c > r d~2 
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where c and d are two constants. Although approaching T c keeping r finite the first 



term is the leading one in the above expression, one finds that the sum in (A.l) is 
dominated by the second term. Within mean field theory (cf) 2 x ) vanishes linearly with 
the temperature at the transition, hence we finds that the disorder affects the critical 
properties for dv u < 2. We verified that this results holds in more general cases like <p n 
field theories and with more general bare propagators. It is natural to conjecture that 
it holds in general above the upper critical dimension. 

In summary, we have found two different Harris criteria. The most important 
physical consequence is that in large enough dimension (d > 2/u u ) critical properties 
of an infinite system will not be affected by an infinitesimal disorder whereas the FSS 
maybe affected even in the infinite dimensional limit depending on the value of the ratio 
d u v u /2. 

In the following appendix we will give a solvable example of the above scenario, 
namely the disordered Blume-Capel model. 

Appendix B. A simple solvable model: the weakly disordered Blume-Capel 
model. 

The relevance of sample to sample fluctuations is quite natural. However, one can be 
surprised that they affect the dynamical finite-size scaling and not the thermodynamics. 
In order to understand in detail, in a concrete example, the general arguments 
formulated above we will consider the following model defined by the Hamiltonian : 

Hbcm = - 2_j ^'S/- S '- s '/ + A 2j &i ~~ h ' 

i,j i i 

SiE {-1,0,1} (B.l) 

the £i = 1 — 5 + \AVs are independent, identically distributed random variables. U are 
i.i.d. random variables normalized in such a way that the second and the fourth moment 
of £i equals one (this is just a convenient but not at all essential choice, see below). When 
the parameter 5 = this is the pure completely connected (or infinite dimensional) 
Blume-Capel model. By increasing 5 one can investigate the role of disorder. It is well- 
known that the pure model, 5 = 0, has a tri-critical point characterized by a correlation 
length exponent v = 1/2 and an upper dimension d u = 3, that is vd u = 3/2 [70J. So 
this is indeed a case where the arguments of the previous section predict that in high 
dimension the critical properties will not be affected by disorder contrary to FSS. Our 
exact solution will confirm explicitly this result. 

The partition function, for a given disorder realization, is given by: 



Sifi) = - '"f 1 + 2e ~" A cosh P' n + ft >f' 



Predictive power of MCT .-Numerical test and Finite size scaling for a mean field spin glass 33 



In order to study the phase transition properties we perform a Landau-like 
expansion of the free energy. One finds at the sixth order: 



i ln(l + 2e-^) + \ 



_ a(/3) ^ 2 



+ 



4!iV 



4 ™ 4 + ^E^ 6 ™ 6 +°(™ 7 ) 



(B.2) 



where 1 — a(ft), 6(/3) and c(/3) are the coefficients of the Landau expansion of the pure 
model. They are given by: 



a(ft) 
6C9) 
c{ft) 



2(3 



-PA 



1 + 2e-P A 



2{3 3 



2(3 5 



-PA _ 4 e -2/3A 

l + 2e-^ A ) 2 

V _ 26e- 2 ^ A + 64e~ 3 ^ A 



(B.3) 
(B.4) 
(B.5) 



(l + 2e-/? A ) 3 

In the pure model, both 1 — a(/3) and b(/3) may vanish and change of sign. The 
cancellation of 1 — a(/3) gives the critical temperature, the sign of b(/3) the order of 
the transition. The results for the pure model are summarized on figure Bl The tri- 
critical point is defined by the simultaneous cancellation of 1 — a{(3) and b(/3), that is 
(A c , T c ) = (2 In 2/3, 1/3). A quick inspection of equation B.2 shows that the pure tri- 




Figure Bl. Phase diagram of the pure fully-connected Blume-Capel model. 



critical point is not suppressed by disorder. Giving a disorder realization, b((3) vanishes 
on the line ft A = In 4: the tri-critical point necessarily sits on this line. More precisely, 
its locus is given by: 

(A c (iV), T C (N)) = (^£^E C 2 ) (B.6) 

i i 

In the thermodynamic limit this gives back the tri-critical point of the pure model: 
T C (N) = T c (oo) = T c . Other choices of normalization of the /j's would have alter 
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the location of the tricritical point but not changed the following results on critical 
properties and FSS. In the thermodynamic limit (large N) ^ Yl^f * s distributed with a 
Gaussian law of width criV -1 / 2 , where a depends on 5. This leads to fluctuations of T c 
of order iV -1 / 2 and, as a consequence, similar fluctuations of A c because of the relation 
(3 A = In 4. This model mimics the ROM case: the thermal fluctuation given by vd u 
are of order N~ 2 ^ 3 (this can be checked both for the pure and the disorder model) the 
disorder fluctuations, given by the fluctuations of the tri-critical point, of order iV" 1 / 2 . 

To study the tri-critical properties of this model, we focus on the line f3A = In 4. 
Writing (3 = f3 c (N) + X/N 2/3 , that is A oc N 2 ' 3 {T C (N) - T), one obtains at leading order 
in N: 

m 




270 

F{x) = Ai 2 (x) + Bi 2 (x) , a = — (B.7) 

8 

where, Ai(x) and Bi(x) are the well-known Airy functions. They are two linearly 
independent solutions of the equation: y" — xy = 0. Their asymptotics, which would 
be useful to obtain the tails of the scaling functions of the different thermodynamic 
observables, is quite simple: 

Ai 2 (x) + Bi 2 (x) — x<-l, 

7TV — X 

± x 3/2 

Ai 2 (x) + Bi 2 (x) ~ 63 1/0 x > 1 . (B.8) 

As a consequence for a given disorder realization, the partition function has a 
standard finite-size scaling form like for equation 21 and that T is independent of the 
realization of the disorder. 



As sketched in subsection |5.2[ averaging on the disorder is here strictly equivalent 
to average on the distribution of tri-critical temperatures: 



For instance, let us detail the computation of the fluctuations of the order parameter 
m defined as: 



m N 



Using the definition of < m >, one gets for a given disorder realization: 
with, 

51840 1 / 3 vr J Bi(«x) + 2^(1; 4/3, 5/3; 15x 3 ) 



/(*) 



a7r 3 / 2 (Ai 2 (xa) + Bi 2 (xa)) 
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And in a more transparent form: 

(T > T c (N))f(x > 1) oc e~ cst * 3/ V /4 

(T < T c (N))f(x < -1) oc (-x) 1 / 4 

As expected, the tail exponent of the scaling function are those one would have find 
if one had computed the magnetization by steepest descent. One see clearly, that the 
magnetization plays in this model the same role as the dynamical overlap. 

Similarly, one can compute the non-connected fluctuations of the order parameter, 

N <™ 2 >= N2,3 my^-MNf ) + 0{l) 



where 

g( x ) 



f ^^yl^Ai' {xa)Ai{xa) + Bi'(xa)Bi(xa) 



Again the asymptotics of g gives back the thermodynamic exponents. 

(T > T c {N))g{x > 1) oc - 

x 

(T < T c (N))g(x < -1) oc ^fx 
The non-connected fluctuations have the same behavior as %4 for the ROM. 

h(x) 




-2-1 12 
Figure B2. Scaling- function of N < Sm 2 >, on this plot x > corresponds to T > T c 



To finish this exercise the connected susceptibility is given by, see figure |B2 
N <5m 2 > = N 2/3 h ' A 



h{\x\ > V 



X 



Now, we have to perform the average over the disorder. For this purpose, let us 
introduce x = y/N Tc ^ N ^~ Tc . One has: 



0(T) 



]\f2a/3 



dxe-^T{N 2/ \e) + N l 



/6 



X) 



(B.ll) 
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As showed in the previous appendix, this integral has two parts: a regular part, 
corresponding to the small value of x, and a singular part given by the tails of T{x) ~ x~ a 
for x » 1. 

To compute this part, one has to consider all the contributions coming from the 
temperature range: \T-T C (N)\ > ciV" 2/3 i.e. such that ax > -N^ 2 (T - T c ) + ciV" 1 / 6 . 
By defining ox Q = N l l 2 {T — T c ), one gets: 

0(T) ~ -= / ^ r ( B - 12 ) 

V27T i CT ( a; +xo)> C 7V-i/6 (x + x) a 



We are now exactly in the same framework as discussed in section 5.2 These simple 
examples show how the analysis of this toy model is instructive. It is a fully solvable 
example of a weak-disordered model, for 5^0, with an exponent v — 1/2 and an 
upper dimension d u = 3, that is vd u = 3/2 [70] where FSS is not the standard one. It 
gives also a clear insight in understanding the competition of the thermal and disorder 
fluctuations. 
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